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Abstract. We use archival Hubble Space Telescope observations of broad-band images from the ultraviolet (F255W- 
^ . filter) through the near infrared (NICMOS F160W-filter) to study the star cluster population of the interacting 

^Nf^ ' spiral galaxy M51. We obtain age, mass, extinction, and effective radius estimates for 1152 star clusters in a region 

, of ~ 7.3 X 8.1 kpc centered on the nucleus and extending into the outer spiral arms. In this paper we present the 

' data set and exploit it to determine the age distribution and relationships among the fundamental parameters 

00 ' (i.e. age, mass, effective radius). We show the critical dependence of the age distribution on the sample selection, 

^«.^ , and confirm that using a constant mass cut-off, above which the sample is complete for the entire age range of 

' interest, is essential. In particular, in this sample we are complete only for masses above 5x10^ Mq for the last 

1 Gyr. Using this dataset we find: i) that the cluster formation rate seems to have had a large increase ~ 50-70 
. Myr ago, which is coincident with the suggested second passage of its companion, NGC 5195, ii) a large number 

CL(' of extremely young (< 10 Myr) star clusters, which we interpret as a population of unbound clusters of which 

JL I a large majority will disrupt within the next ~10 Myr, and Hi) that the distribution of cluster sizes can be well 

^ ■ approximated by a power-law with exponent, — r; — —2.2 ± 0.2, which is very similar to that of Galactic globular 

clusters, indicating that cluster disruption is largely independent of cluster radius. In addition, we have used this 
, dataset to search for correlations among the derived parameters. In particular, we do not find any strong trends 

between the age and mass, mass and effective radius, nor between the galactocentric distance and effective radius. 
There is, however, a strong correlation between the age of a cluster and its extinction, with younger clusters being 
, more heavily reddened than older clusters. 

^ '■ 
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1. Introduction 

This study aims at understanding the formation history 
of star clusters, their properties and spatial distribution in 
the interacting spiral galaxy M51. To study such a system, 
one needs the superb spatial resolution of the Hubble Space 
Telescope (HST) in order to avoid crowding effects and 
to differentiate between individual stars, associations, and 
compact star clusters. 

Although much work has already been done on extra- 
galactic star clusters using the HST, most studies have 
concentrated on specific components of the full cluster 
populations (e.g. star clusters in the center of spiral galax- 
ies, Boker et al. 2001) or star cluster systems in extreme 
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environments such as galactic mergers (e.g. Miller et al. 
1997, Whitmore et al. 1999). In order to form a baseline to 
study the effects of environment on cluster formation, evo- 
lution and general cluster properties, one must study the 
properties of clusters in normal (i.e. more common) envi- 
ronments. It is only then that we are able to see which 
properties are truly unique for a given environment and 
which properties remain fixed across all environments. 

One such property of star cluster systems that has a 
well-established baseline is that of the luminosity function. 
Larsen (2002) has shown that the cluster populations of 
spiral galaxies follow a luminosity distribution that is well 
represented by a power-law, N{L) oc i^", with q:=2.0. 
This is remarkably similar to that found in all other en- 
vironments (de Grijs et al. 2003c), and hence can be re- 
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garded as a general property of young star cluster systems. 



Recent work by Larsen (2004) shows that young clusters 
in spiral galaxies have typical sizes of R^s ~ 3-10 pc, with 
the youngest clusters having extended halos. It remains 
to be seen how other properties, such as mass, age and 
spatial distributions depend on the environment in which 
the clusters form. 

M51 provides an almost ideal astrophysical laboratory 
to study extragalactic star clusters. This is due to its rela- 
tively close distance of ^ 8.4 Mpc (Feldmeier et al. 1997) 
and its almost face-on orientation. Physically it is an in- 
teresting case study because it seems to have had a strong 
interaction with its companion, NGC 5195 an SO pecu- 
liar galaxy, during the last few hundred Myr, which pre- 
sumably caused its grand design spiral appearance as well 
as its high star formation rate. The M51/NGC 5195 en- 
counter(s) have been modelled in great detail by many 
authors; this will allow us to compare the derived age dis- 
tribution of the star clusters with the orbital parameters 
of the system, which are robustly established. The system 
has been modelled relatively successfully by a single early 
passage --300 ± 100 Myr ago (Toomre & Toomre 1972, 
Hernquist 1990, and Salo & Laurikainen 2000, hereafter 
SLOO). Additionally, SLOO propose a model of a double 
passage, the first happening ~400-500 Myr ago and the 
second, or last encounter, happening ~50-100 Myr ago. 
This latter model seems to fit many observed details that 
the single passage models fail to reproduce, such as the 
observed counter-rotation of the southern Hi tail. 

Much work has already been done on the star cluster 
system of M51. A full literature review is beyond the scope 
of this paper, but we refer the reader to the following work: 
Bik et al. (2003, hereafter Paper I) have studied the star 
cluster population in a relatively small area to the north 
east of the nucleus that included inner sections of one of 
the spiral arms. They find that the mass function is rea- 
sonably well fit by a power-law of the form N{M) oc M~" 
where a ^ 2.0. They also show that there is weak evi- 
dence for a possible increase in the cluster formation rate 
in the inner spiral arms at ~ 400 Myr ago, the proposed 
time of an interaction between M51 and NGC 5195, but 
that the evidence is only at the 2 a level. Boutloukos & 
Lamers (2003) used the data from Paper I to determine 
the characteristic disruption timescale of clusters in this 
region, which is ^ 40 Myr for a \Q^Mq cluster. Larsen 
(2000) has studied a few of the "young massive clusters" 
which form the high end of the luminosity function, using 
combined HST and ground-based observations. He finds 
that although M51 has a relatively high star- formation 
rate, it fits the specific L/-band luminosity (T^) vs. star 
formation rate relation of non-starburst galaxies. This has 
been interpreted as showing that high-mass star clusters 
will form whenever the star formation rate is high enough. 
Scoville et al. (2001) have used many of the same HST 
images used in this study, to look at OB star formation 
in M51. They find that the mass function of OB star 
clusters is well represented by a power law of the form 



^(Mci)/dMei cx 



-2.01 



in agreement with the results 



presented in Paper I. They also show that although the 
spiral arms only make up ~ 25% of the disk surface area 
they contain a much higher fraction of the catalogued Hii 
regions. 

This study is the second in a series of three, which aim 
at understanding the properties of the star cluster system 
of M51 as a whole. This paper introduces the data and 
techniques used to determine the cluster properties, anal- 
yses the age distribution as well as the correlations among 
the cluster parameters, namely their age, mass, extinction, 
size, galactocentric distance, and position in the galaxy. In 
the third paper IjGieles et al. 20041 hereafter Paper HI) we 
determine the luminosity and mass distributions and de- 
rive the disruption timescale for clusters in M51. 

The study is setup in the following way: in § |21 we 
present the observations used in this study, including the 
reduction and photometry. In § |3 the methods of find- 
ing the clusters and the selections applied to the data are 
elucidated. In § 01 we present our method of determining 
the ages, extinctions, and masses of each of the detected 
clusters, and discuss how contaminating sources are de- 
tected and removed. In §§ Eland El the general properties 
of the cluster system are presented and discussed, while 
§13 presents evidence for a young, short lived cluster pop- 
ulation. In §§ [HI and ini the size distribution and relations 
between size and the other cluster parameters are pre- 
sented, followed by 1101 in which we summarize our main 
results. 



2. Observations, reduction and photometry 

Our analysis of stellar clusters in M51 is based on ob- 
servations done with the Wide Field Planetary Camera- 
2 (WFPC2) and the Near-Infrared Camera and Multi- 
Object Spectrometer (NICMOS) on board the HST. 



2.1. WFPC2 Observations 

Observations in 8 different passbands of the inner region 
of M51 were obtained from the HST datB, archive: F255W 
(« UV), F336W (« C/), F439W (« B), F502N (w OIH), 
F555W (« V), F656N (« Ha), F675W (w R) and F814W 
(« /). We emphasize that we have not transformed the 
HST filters into the standard Cousins- Johnson filter sys- 
tem. 

An overview of the filters and exposure times is given 
in Table 1. Proposal ID 5777 and 7375 contain the F555W 
band images and will be referred to as Orientation 1 and 
Orientation 2, respectively. All the other orientations are 
rotated/shifted to these 2 base orientations. An overview 
of the orientation of the data sets is shown in Figures ^ 
and 121 The combined data sets span a region of about 
180" X 200" (i.e. 7.3 x 8.1 kpc) containing the complete 
inner spiral arms of the galaxy. 

The raw images were flat fielded using the auto- 
matic standard pipeline reduction at the Space Telescope 
Science Institute (STScI). Warm pixels were removed us- 
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ing the warmpix task which is part of the IRAF/STSDAS^ 
package. Bad pixels where fixed with the wfixup task 
which interpolates over bad pixels in the x-direction. 

Cosmic rays where treated differently for the three 
datasets. For the images consisting of multiple exposures 
of the same field per filter the task crre j was used. This 
combines exposures of the same field by rejecting very high 
counts on an individual pixel basis. Three iterations were 
done with rejection levels set to 8, 6 and 4(t. 

In single exposure images the cosmic rays were re- 
moved using the Lacos_im package of van Dokkum (2001) 



This algorithm identifies cosmic rays of arbitrary shapes 
and sizes by the sharpness of their edges and reliably dis- 
criminates between poorly sampled point sources and cos- 
mic rays. The parameter setting used was a sigclip of 
6.5, a sigf rac of 0.5 and an objlim of 4. Four iterations 
were done to give the best result, because the result of 
the cosmic ray rejection is very sensitive to the settings. 
Careful tests were done to make sure no real sources were 
removed by the cosmic ray removal task. This resulted in 
almost cosmic ray free images. The few remaining cosmic 
rays will be rejected later when the photometry coordi- 
nates from the source detection in the different filters will 
be cross correlated. 

2.2. NICMOS Observations 

NICMOS images were obtained with the NIC3 camera, 
which has a field of view (FoV) of about 52"x52". A mo- 
saic of nine overlapping pointings of the NIC3 camera of 
M51 was available in the archive, yielding a total FoV of 
186"xl88" (i.e. 7.6 x 7.7 kpc) of the central region in the 
near-infrared FllOW and F160W filters (see Fig.[2). Each 
of the nine mosaic positions was observed using a square 
dither in each filter setting. 

Flat-fielding, dark frame correction and bias subtrac- 
tion were all performed in the pipeline image reduction 
and calibration by the calnica task. Bad pixels were re- 
paired with the f ixpix task which uses static bad pixel 
mask images. All individual images in the dither pattern 
were shifted to the same orientation with imshift. The 
dither pattern constituted of shifts of an integer number of 
pixels, the images were perfectly aligned on subpixel level. 
The individual images in the same orientation were added 
with imcombine using the crrej method for cosmic ray 
rejection. Finally, the nine co-added images were mosaiced 
by matching sources in the overlap region and comparing 
histograms of the sky values. This results in a large square 
FoV overlapping a large area of the optical data set. Since 
the mosaic routine does not scale the count rates values 



^ The Image Reduction and Analysis Package (IRAF) is dis- 
tributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. STSDAS, the Space 
Telescope Science Data Analysis System, contains tasks com- 
plementary to the existing IRAF tasks. In this study Version 
2.1.1 (December 1999) was used. 



when matching the individual chips, the mosaic images 
can be used for source selection and photometry. 

2.3. Source selection 

Point-like sources were identified in the individual WFPC2 
chip images and the NIC3 mosaic images with the daof ind 
task from the DAOPHOT package Stetson (1987). For the 
optical filters a threshold of 5a was used, where a is the 
sky noise. In order not to be biased towards certain ar- 
eas on the chips, the value of a was determined on a low 
background region on each of the chips. For the F255W 
filter the threshold was Aa and for the near-infrared NIC3 
images the threshold was set to 8a. The threshold was cho- 
sen such that in each filter and chip the resulting source 
density was more or less the same. No limitation was set 
on sharpness or roundness, to include as many sources as 
possible, including unrelaxed young clusters. 

2.4. Photometry 
2.4.1. WFPC2 

The coordinates from the source lists, obtained from the 
daof ind output, were used as an initial guess of the cen- 
ters for the phot aperture photometry routine. For the 
WF chips we used an aperture, annulus (inner boundary 
of background annulus) and dannulus (width of the back- 
ground annulus) of 3, 7 and 3 pixels respectively. For the 
PC chip we used slightly larger values, namely 4, 10, 3 for 
the aperture, annulus and dannulus respectively. The cho- 




13''30"'00= 13''29"'45' 13* 

Right Ascension (J2000) 

Fig. 1. The orientation of the two iJST pointings for the 
F336W, F439W, F555W, F675W, F814W, and F656N fil- 
ters overlaid on a Digital Sky Survey image of M51. North 
is up, and east is to the left. The left orientation/pointing 
is referred to as field 1, while the right is referred to as 
field 2. 
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Table 1. Overview of the datasets used. ^Equivalent to the V3 Position Angle of the WFPC2. 



ID 


Filters 


Exp. time 


Observation Date 


Orientation angle^ 


7375 


F336W 


2 X 600 s. 


21 July 1999 


275.9 


(WFPC2) 


F439W 


500 + 600 s. 








F555W 


2 X 600 s. 








F656N 


1300 + 700 s. 








F675W 


500 s. 








F814W 


700 + 300 s. 






5777 


F439W 


2 X 700 s. 


15 Jan 1995 


101.5 


(WFPC2) 


F555W 


600 s. 








F675W 


600 s. 








F814W 


600 s. 






5123 


F502N 


400 + 1400 s. 


24 January 1995 


93.2 


(WFPC2) 


F656N 


500 + 1200 s. 






5652 


F255W 


4 X 500 s. 


12 May 1994 


333.0 


(WFPC2) 


F336W 


2 X 400 s. 






7237 


FUOW 


128 s. 


28 June 1998 


-113.2 


(NICM0S3) 


F160W 


128 s. 







sen values for the apertures are ~ 3 times larger than the 
FWHM of the PSF in order to account for the fact that 
at the distance of M51, clusters appear slightly extended 
(i.e. are not point sources). Photometric calibration was 
done by applying zero-point offsets from Table 28.2 of 
Volt (1997)1 in the VEGAMAG system. 

Total magnitudes were determined by correcting for 
flux outside the measurement aperture. Ideally aperture 
corrections are determined from real sources on the im- 
age. In the case of our data there were not enough iso- 
lated bright sources located in regions of low background 
to obtain a meaningful correction. Thus aperture correc- 
tions from the measurement aperture radius to a 0.5" ra- 
dius were determined from analytically generated clus- 
ters by convolving a PSF (generated using Tiny Tim, 
Krist & Hook 1997) with King profiles using the Baolab 
package (Larsen 1999). King profiles of concentration fac- 
tor, c, of 30 (c = rt/rc where rt is the tidal radius and 
Tc is the core radius), and effective radius, i?cff, of 3 
pc (see § IHJl were used for the analytic cluster profiles. 
The values for the three WF chips are averages. For the 
PC and WF chip in each filter the values are listed in 
Table 3. A final —0.1 mag. correction was applied to all 
sources to correct for the light missed in the 0.5" aperture. 
After applying the aperture corrections, the magnitudes 
where corrected for CTE loss according to the equations 



of Whitmore et al. (1999a) 



2.4.2. NICMOS 

The photometry procedure for the NIC3 images is essen- 
tially the same as for the WFPC2 images. An aperture, 
annulus and dannulus of 2, 5 and 3 pixels were used, re- 
spectively. The fiux was converted to magnitudes in the 
VEGAMAG system using the PHOTFNU values from 
the headers. Aperture corrections from the measurement 
aperture radius to a 1.0" radius were determined using 



the same technique as for the WFPC2 images. Values are 
listed in Table 2. A final correction to a nominal infinite 
aperture correction of —0.08 mag (i.e., multiplying the flux 
with 1.075) was applied. 

2.5. Incompleteness due to detection limits 

In order to quantify at what magnitude our sample starts 
to be incomplete due to detection limitations, the 90% 
completeness limit was determined for each band for each 
pointing. To this end, artiflcial sources were added to the 
HST images and were recovered with exactly the same 
routines as used to find the real sources. A very impor- 
tant criterion is that objects should be 5 times above the 
standard deviation of the background. The value for a 
should be the same as when real sources are found. The 
magnitude of the added sources were determined and an 
error criterion of Amag < 0.2 was applied (see §|S1), the 
same criterion as for the real sources to enter the photom- 
etry list. The tests were done for the PCI and WF3 chips 
of both pointings separately. 

Extended sources with different i?cff are simulated by 
convolving WFPC2 PSFs, produced by Tiny Tim, with 
King profiles. The King profiles are produced by the 
Baolab package (Larsen 1999). King profiles with a con- 
centration index, c, of 30 with different i?eff were tried. In 
§|Hlit is shown that the average size for resolved clusters in 
M51 is ~ 3 pc. Thus, we take a more conservative estimate 
of the completeness survey by assuming i?cff = 5 pc (the 
more extended the cluster is the more diffuse it becomes, 
for a given brightness, making it more difficult to detect). 
On each chip, two regions were studied, one with low and 
one with high background levels. Since the aperture cor- 
rections were derived for average extended sources (i?off — 
3 pc), the magnitudes of the artificial 5 pc sources are too 
faint and the magnitudes of the pure PSFs are too bright. 
In both cases the maximum offset is about 0.15 mag. 
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The recovered fraction as a function of magnitude in 
the F555W filter is plotted in Fig. for the WF3 chip. 
Each plot shows the effects of background and source size. 
The plots show that for an individual chip the 90% com- 
pleteness limit varies by more than a magnitude, depend- 
ing on the location on the chip and the size of the source. 
We therefore do not have a clear detection limit but a re- 
gion where incompleteness sets in. Details of this effect and 
its influence on the luminosity function are given in Paper 
III. The worst case 90% completeness limits are shown in 
Fig.0]for all bands and for the 3 available ifS'T pointings. 
When more than one pointing is available for a filter, the 
brightest (i.e. most conservative) 90% completeness limit 
of that filter is chosen. 

For the remainder of this study the worst case 90% 
completeness limit was applied to the data. Since the PC 
chip has been shown to have many single bright stars 
(Lamers et al. 2002) and since its completeness limit is 
about 1 magnitude brighter than for the WF3 chip, we 
will only use the WF data for the remainder of this study. 

3. The selection of the clusters and the 
photometry master list. 

The coordinate lists for the individual WFPC2 chips were 
transformed to one system using the metric task. This 
task corrects for geometric distortion and yields coordi- 
nates with respect to one frame (WF3). Transformation 
matrices to convert the mosaic coordinates of field id 
#5652 (F255W and F336W) (see Table 1) to Orientation 
1 and Orientation 2 coordinates, were determined by find- 
ing sources in common to both frames in the F336W fil- 
ter. Since the #5777 program does not contain an F336W 
exposure, the F439W exposure was used instead. The 
sources used to find the transformations were selected to 




Reft = 5 PC background = high 

Reft = 5 pc , bockqround = low 

Reff = 0, background = high 

Re(f = 0, background = low 



Fig. 3. Recovered fraction of sources added to the F555W 
science images on the WF3 chip. Four examples are shown, 
combinations of low and high background and PSFs and 
extended clusters {ReS = 5 pc). 



WF3 



+ Field 1 
O Field 2 



E 

E 20 - 



F439W F555W 
Filter 




Right Ascension (J2000) 



Fig. 2. The orientation of the F255W FoV (black) and the 
NIC3 mosaic (white). 



Fig. 4. Worst case (high background, Rcs = 5 pc clusters) 
90% completeness limits for PCI and WF3 for all filters 
and the three available HST pointings. 



be isolated and spread over the images by as much as 
possible to make the transformation as accurate as possi- 
ble. The spatial transformation functions were calculated 
based on the coordinates of the sources in the two frames 
using the geomap task in IRAF. The accuracy of the trans- 
formation is approximately 0.1-0.2 pixels. The same was 
done for the near-infrared filters (FllOW and F160W). 
For these transformations identical sources were found in 
the FllOW and F814W bands. In this case the accuracy 
was somewhat lower since the NIC3 pixel size (wO.2") is 
larger then the WF pixel size (wO.l"). The accuracy of 
the transformation function is still 0.2-0.3 WF pixels. 

The final source list was made by cross correlating 
sources found in the various filters. First, a cross corre- 
lation between the F439W, F555W and F675W was done. 
If a coordinate was found in all 3 filters, allowing a mis- 
match of 1.4 pixels, it was defined as a genuine source. 
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Table 2. 90% completeness limits for the WF3 in the 3 
different available pointings. Most conservative values are 
printed bold. 



PCI WF3 



Filter 


1 


2 


3 


1 


2 


F255W 






18.80 






F336W 




20.8 


20.45 




21.60 


F439W 


22.10 


21.80 




22.90 


22.60 


F555W 


21.80 


21.90 




22.70 


23.20 


F675W 


21.20 


21.90 




22.40 


22.45 


F814W 


21.20 


21.00 




22.20 


22.20 



adopting only foreground extinction of Ay — 0.40 mag 
and E{B - V') = 0.13 mag (Massey et al. 1995) and a 
distance modulus of 18.48 mag (Westerlund 1997) The 
total number of stars brighter than My < —7.0 mag 
(which corresponds to V < 22.60 at the distance of M51, 
for no extinction) is less than 100, and the number with 
Afy < —8.0 mag is 6. By far the majority of these stars are 
found in clusters. We can get an estimate of the fraction 
of the luminous field stars compared to the total fraction 
from the data in Table 12 of Massey (2002) which gives the 
number of massive stars outside clusters. (This is the only 
large and homogeneous extragalactic stellar sample deep 
enough to determine the relation between the populations 
of field stars and the total populations.) He found 186 field 
stars (i.e. outside known OB associations) more massive 
than 40 Mq, which corresponds to Mboi < —9 mag. More 
than 90 per cent of these are on or near the main sequence, 
where Tcs > 25 000 K and hence the bolometric correction 
is so large that My > —7.0 mag. At most 10% of these 
stars, i.e. less than 18, are optically bright at My < —7.0 
mag in the LMC. 

With a total mass of the observed area of about 4 x 
10^ Mq (Meatheringham 1991), this implies about 4 to 5 
field stars more massive than 40 Mq per 10^ Mq. (This 
number could be slightly higher if the dark matter content 
of the LMC within the region used here is considerable). 
This small number is due to the very steep IMF of the 
massive field stars. Massey (2002) showed that the stellar 
IMF of massive field stars can be approximated with F ~ 
—4, compared to F ~ —1.3 for cluster stars, where F = 
l + d(log7V(Af))/d(log(M)). 

The total luminous mass of M51 within a radius of 5 
kpc is 2.1 X lO^^M© (Salo & Laurikainen 2000). Scaling 
the number of field stars per unit mass in the LMC to 
that of M51, we estimate that at most 170 to 210 field 
stars with My < —7.0 will be in our field, and less than 
15 to 20 with My < —8.0, if there were no extinction in 
M51. Even a moderate extinction will remove these stars 
from our sample. 

Another estimate of contaminating stellar sources in 
our sample can be obtained by comparing the number of 
field stars to that of observed clusters. This comparison 
relies on the assumption that the ratio of field stars to 
clusters is the same in different galaxies. This assump- 
tion will not be valid for low or intermediate mass stars 
between galaxies with widely different cluster disruption 
timescales, as galaxies with shorter disruption timescales 
will have a higher field star to cluster ratio. However, for 
young massive stars (lifetimes less than 10 Myr), disrup- 
tion is not expected to significantly influence this ratio. 
Using the catalog of Bica et al. (1996) we see that there are 
~ 180 clusters and associations brighter than My = —6.94 
in the LMC (assuming a distance modulus to the LMC 
of 18.48 and an imposed completeness of my — 22.70 at 
the distance of M51). Using the arguments given above, 
there are < 18 young massive stars that are brighter than 
My — —7.0 in the LMC. Thus there are less than one field 
star for every ten clusters in the LMC. Therefore, using 



Table 3. Aperture correction values for all filters (mag). 



Filter 


PC 


WF3 (mean) 


NIC3 


F225W 


-0.32 


-0.14 




F336W 


-0.31 


-0.13 




F439W 


-0.31 


-0.13 




F555W 


-0.32 


-0.13 




F675W 


-0.32 


-0.13 




F814W 


-0.33 


-0.14 




FllOW 
F160W 






-0.22 ±0.01 
-0.22 ±0.01 



' This value is applied to our sample. For clusters with Rcb = 5 
pc, these values decrease by ^ 0.15 mag. (e.g. to -0.28 for 
F555W in the WF3 chip). 

Then the other filters were compared with the source list, 
again allowing a mismatch of 1.4 pixels. The transformed 
coordinates of the F255W and F336W frame were allowed 
to have a 1.6 pixel mismatch. The FllOW and F160W 
frames were allowed to have a 1.8 pixel mismatch. This 
was adopted because of the uncertainty in the transfor- 
mation. 

The resulting initial source list contains 3504 objects 
that are observed in at least the F439W, F555W, and 
F675W bands above the 5a detection limit. Most of these 
will be eliminated later in the study of their energy dis- 
tributions because the photometry is uncertain or the en- 
ergy distribution does not match that of cluster models 
(see 

3.1. Contamination of the sample by stars 

We estimate the possible number of massive stars that 
may contaminate our source sample by scaling the star 
population of the LMC, measured by Massey (2002) to 
the conditions in M51. This sample, which is complete 
down to y ~ 15.7 mag or My ~ —3.2 mag covers a to- 
tal area of 10.9 kpc^ and includes most of the bar, the 
30 Doradus region and several other well populated fields. 
This area contains more than half of the mass of the LMC. 
We have calculated the intrinsic visual magnitudes, by 
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this approximation, we expect less than 10% of our M51 
sample to be bright young stars. 

A third alternative method to estimate the contamina- 
tion by massive isolated field stars is to scale the number 
of isolated field stars in the LMC to that of M51 using 
the mass normalized star formation rates. In the studies 
summarized by Grimm et al. (2003), M51 has about half 
the star formation rate per unit galaxy mass as the LMC. 
So, using this assumption we expect half as many bright 
field stars in M51 as in the LMC per unit mass, which 
results (using the mass ratios above) in ~ 47 bright field 
stars. 

We conclude that the number of contaminating stars 
in our sample is small^. Moreover, we will show in Section 
14.51 that most of these will be removed from the sample 
because they do not pass the criterion that their energy 
distribution can be fitted with a cluster model. 

4. Determination of the cluster parameters from 
their energy distribution 

The age, mass and extinction of star clusters can be de- 
rived by comparing their spectral energy distributions to 
those of cluster evolution models (e.g. Paper I, de Grijs et 
al. 2003a, b,c). This method can also be used to eliminate 
contaminating sources (e.g. stars and background galax- 
ies). 

4.1. Adopted cluster models 

We have adopted the updated GALEV simple stellar pop- 
ulation (SSP) models (Schulz et al. 2002, Anders & Fritze 
v. Alvensleben 2003) as our templates to derive ages, ex- 
tinction values, and masses for each of the observed clus- 
ters in M51. These models use the isochrones from the 
Padova group, which seem to fit the data better than the 
Geneva tracks (Whitmore & Zhang 2002). The set of mod- 
els that we use assumes a Salpeter stellar IMF from 0.15 
- 50 Mq. We have used models with metallicities of 0.4, 
1, and 2.5 Zq, and a comparison between the different 
metallicities will be made in Section ITHl 

There were two major reasons for selecting the up- 
dated GALEV tracks. The first is the inclusion of gaseous 
emission lines as well as continuum emission, which can 
severely affect the broad-band colours for young star 
clusters (see Anders & Fritze-v. Alvensleben 2003). The 
second is that they have been published in the HST 
WFPC2/NICM0S filter system, which eliminates possi- 
ble inaccuracies in converting the _ff5'T observations to the 
standard Cousins- Johnson system. 

^ Bik et al. (2003) found that massive young clusters in the 
inner NE-spiral arm of M51 have a surprisingly small flux in 
the 0[III] A5007 line. They interpret this as evidence that there 
is a lack of massive stars of M > 30Mq in the clusters. If this 
is also valid for the field stars, the number of contaminating 
field stars in our sample will be even smaller than estimated 
here. 



4.2. The fitting of the spectral energy distributions 

We used the three-dimensional maximum likelihood fit- 
ting method {3DEF) developed by Bik et al. (2003). The 
method has been described in detail there, and hence we 
will only give a brief summary of it here. In passing we note 
that this method has been tested against colour-colour 
methods of age fitting, and has shown itself to be supe- 
rior (e.g., de Grijs et al. 2003a vs. Parmenticr, de Grijs & 
Gilmorc 2003). 

The method works as follows: we use a grid of SSP 
models, i.e. the GALEV models (see previous section) that 
give the broad-band colours of a single age stellar popu- 
lation as a function of age. For each model age, we apply 
50 different extinctions in steps of 0.02 in E{B — V). We 
adopted the Galactic extinction law of Savage and Mathis 
(1979), which was found to agree with the extinction law 
derived for M51 by Lamers et al. (2002). We then com- 
pare the model grid with the observed spectral energy 
distribution of the clusters using a minimum test, with 
each observation weighted with its uncertainty. The model 
(age and extinction) with the lowest is selected, and 
the range (i.e. maximum and minimum age) of accepted 
values is calculated by taking the most extreme model 
that satisfies xt < xtmin + 1- Once the age and extinc- 
tion have been calculated, the mass is found by comparing 
each observed filter magnitude to each model magnitude 
(iteratively with respect to age and extinction). Since we 
have first scaled the GALEV models to a present mass of 
1O^M0 the difference between the observed and the pre- 
dicted magnitude (taking into account the distance mod- 
ulus and extinction) can be converted to the present mass 
of the observed cluster^. 

4.3. Test of the accuracy of the method 

We have tested the accuracy of the method used for the 
determination of the age and mass of clusters. To this 
purpose we created a grid of synthetic clusters, distributed 
evenly in logarithmic age, in the range of 7.0 < log(t/yr) < 
10.0 in 31 steps, and evenly in logarithmic mass in the 
range of 3.0 < log(M/Mo) < 6.0 in 21 steps. This re- 
sults in 651 synthetic clusters. The clusters have a dis- 
tance modulus of m ~ M = 29.64, as for M51, and were 
assigned a random extinction of < E(i? — V) < 1.0 
mag with a Gaussian probability function that peaks at 
E(i? — V) — mag and has a a = 0.20 mag. For each 
of these clusters we calculated the expected energy distri- 
bution in the same HST/WFPC2 filters as used for this 
study of M51, using the GALEV cluster models with a 
metalhcity of Z = 0.02, and a Salpeter IMF. The pho- 
tometry of each synthetic cluster was given an uncertainty 

^ AU masses derived in this paper refer to the present mass 
of the clusters, not to their initial mass. If stellar evolution is 
the only mass loss mechanism of the cluster, the initial mass 
of the clusters can be retrieved by correcting the present mass 
for the fraction lost by stellar evolution, given by the GALEV 
models. 
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Table 4. Empirical parameters for the uncertainty of the 
magnitudes: A(m) = lO'^i+'^^xm 



Filter 


di 


d2 


F255W 


-8.540 


+0.390 


F336W 


-7.850 


+0.325 


F439W 


-8.320 


+0.328 


F555W 


-7.800 


+0.300 


F675W 


-7.600 


+0.300 


F814W 


-8.230 


+0.328 



that depends on magnitude. The uncertainty was derived 
from the observations, which show that A(m) can be ap- 
proximated quite well by A(m) = iQd.i+d2xm ^ where m is 
the magnitude and the values of di and o?2 are empirically 
determined for each HST filter. These values are listed in 
Table 01 We then added a random correction to the mag- 
nitude of each cluster in each filter, within the interval 
(— A(77i), +A(m)). This resulted in a list of 651 synthetic 
clusters with photometry and photometric errors similar 
to those of the observed clusters in M51. 

We then ran these synthetic clusters through our clus- 
ter fitting program to determine the age, mass and extinc- 
tion in exactly the same way as done for the observed clus- 
ters. We adopted the same conservative magnitude limits 
as for the observed sample, as given in Tabled The com- 
parison between the input mass, age and extinction of the 
synthetic clusters and the result of the fitting gives an in- 
dication of the accuracy and reliability of the results. The 
result is shown in Figure |31 

The panels show the difference between input and out- 
put values of log(t), log(A/) and E(i? — V) as a. function of 
the V magnitude of the synthetic cluster. This parameter 
was adopted because the accuracy of the fitting is expected 
to depend on the accuracy of the input photometry, for 
which is a good indicator. We have divided the results 
into four age bins (with parameters indicated in the fig- 
ure), because the accuracy of the fit may be different in 
different age intervals, as the spectral energy distribution 
depends on age. 

Figure shows that the age determination is quite 
accurate, to within about 0.20 dex, for relatively bright 
clusters with V < 21.0 mag. For fainter clusters with 

V > 21.0 mag the uncertainty increases up to a maxi- 
mum of about 1 dex for clusters with V > 22.0 mag. The 
scatter is approximately symmetric around the zero-line 
with two exceptions. 

(i) For the youngest faint clusters with i < 30 Myr and 

V > 21.5 mag the age is more likely to be overestimated 
than underestimated. This is the consequence of the sim- 
ple fact that the two youngest GALEV cluster models 
have ages of 4 and 8 Myr. So a cluster cannot be fitted 
to a model younger than 4 Myr, which limits the possible 
negative values of log(iout/iin)- 

(ii) Faint clusters with ages of log(i/yr) > 7.5 have a 



higher probability to be fitted with a younger cluster 
model than with an older cluster model. This might re- 
duce the number of old clusters derived from the fitting 
of the SEDs of observed clusters. This effect will be taken 
into account in the analysis of the observed mass and age 
distributions of the M51 clusters. 

Figure [Sja shows that the mass determination is quite 
accurate for relatively bright clusters with V < 21.0 mag. 
For fainter clusters with V > 21.0 mag the uncertainty 
decreases up to a maximum of about —0.5 dex to +1.0 
dex for the faintest clusters. The scatter is approximately 
symmetric around the zero-line, except for the tendency 
to overestimate the mass of the faint old clusters. This is 
a direct result of the underestimate of the age. 

Figure \^ shows that the extinction can be derived 
quite accurately for young clusters with log(t/yr) < 8.5 
and V < 21.5 mag. For fainter clusters there is a symmet- 
ric scatter up to about AF,{B — V) ~ 0.1 mag. For clusters 
older than 300 Myr there is a tendency to overestimate the 
extinction. This is the reason for the overestimate of the 
mass and the underestimate of the ages of these old clus- 
ters. 

In general the ages and masses derived from fitting 
the cluster models are quite accurate: more than 99% of 
our synthetic clusters have an age determination that is 
accurate to within 0.25 dex and a mass determination that 
is accurate to within 0.2 dex. For real cluster samples the 
accuracy will depend on the magnitude distribution of the 
clusters, and hence on their photometric accuracy. 

4.3.1. Stochastic colour fluctuations 

An additional uncertainty in the fitting of ages of star 
clusters using SSP models, is the stochastic colour fiuctu- 
ations resulting from the fact that young low mass clus- 
ters will be dominated by a few high luminosity stars (e.g. 
Dolphin & Kennicutt 2002). Thus, the absence (or pres- 
ence) of a small number of stars can significantly affect the 
observed integrated colours of a cluster. Unfortunately, 
due to the random nature of this effect it is impossible to 
determine which clusters in our sample are most heavily 
affected by this. However, if the stochastic colour fiuctu- 
ations vary equally around the 'true' (i.e. fully sampled, 
high mass limit) colour, we do not expect our derived age 
distribution for the population as a whole to be affected 
by this. Finally, we note that our age distribution analy- 
sis (see § has been restricted to clusters with derived 
masses above 10^' ''M©, e.g. the most massive clusters, 
where the stochastic colour fluctuations are expected to 
be smallest (Langon & Mouhcine 2000). Stochastic effects 
are expected to affect the colour and magnitude of a 10 
Myr, IO'^Mq cluster by ay = 0.08 mag and ay^i = 0.07 
(Dolphin & Kennicutt 2002). These small fluctuations are 
within the observational uncertainties, and as such are not 
expected to signiflcantly affect the derived parameters. 
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Fig. 5. Comparison between input and output values of 
the ages (upper panel, a), masses (middle panel, b) and 
extinction values (lower panel, c) of the synthetic cluster 
sample. Different symbols indicate different age ranges. 



4,3,2, Dependence on the adopted models 

The 3DEF method assumes that the adopted models accu- 
rately reproduce the colours of star clusters as a function 
of age. In order to quantify this effect we have carried out 
the same analysis as discussed in § 14.31 but now using the 
Starburst99 models (Leitherer et al. 1999) for the artificial 
clusters, and then carrying out our fittting procedure using 
the GALEV models. These models were chosen because 
of the differing stellar evolutionary isochrones used when 



constructing their respective SSPs (i.e. the GALEV mod- 
els adopt the Padova isochrones, while the Starburst99 
models adopt the Geneva isochrones). 

These tests resulted in 80% of the artificial clusters be- 
ing fit within 0.5 dex of the input age and mass. The ex- 
tinction is fit to within 0.2 mag in E(B-V) for > 80% of the 
artificial clusters. Thus we see that while there are signifi- 
cant differences between the models used, the parameters 
of the majority of clusters can be accurately recovered. 

We also note that fitting our cluster sample using the 
Starburst99 models as templates, resulted in age and mass 
distributions which are consistent with the results pre- 
sented here (Bastian & Lamers 2003). Namely, the high 
numbers of very young (< 10 Myr) as well as the increase 
the cluster formation rate at ^ 50—100 Myr ago (see §|H1). 

4.4. Reduction of the sample based on energy 
distributions 

Once we have obtained the best fit for each source (i.e. 
age, extinction and mass) we can then determine whether 
or not the object has a good fit, i.e. if it is a cluster. In 
previous studies (Paper I, de Grijs et al. 2003a) we used 
the selection criterion that the reduced must be below 
3.0. This was effective in eliminating most of the poor fits, 
but it also removed many of the brightest clusters because 
their photometric error was relatively low. In order to re- 
move this problem, we have used the standard deviation, 
defined to be 



a\BVR) 



E 



(m: 



mod 



nmte 



(1) 



where the sum runs over the F439W, F555W, and F675W 
bands only (the three bands in which all clusters were 
detected), as a criterion to test the accuracy of the photo- 
metric fit. The sum was restricted in order to maintain a 
fairly homogeneous set and to eliminate biases that would 
occur if bands with short exposures or low count rates 
(e.g. F255W) were included. 

We then checked by eye the spectral energy distribu- 
tions of the sources and their fits, in order to find a reason- 
able value for the minimum a^{BVR), o'min(SV^^)- As the 
"goodness" of the fit is rather arbitrary, we quantified the 
exact value of a'^i^{BVR) by making a histogram of all 
the sources, approximating the histogram with a Gaussian 
(which is a relatively good fit) and taking the 1/e point. 
This point, and the value preferred by our visual inspec- 
tion were almost identical at a^^^^^{BV R) — 0.05. 

In the remainder of the paper, when we refer to 
"clusters" , we will mean sources that passed the following 
criteria: 



detected in the F439W, F555W, and F675W bands at 

at least 5ct above the background, 

photometric uncertainty < 0.2 mag in the F439W, 

F555W, and F675W bands, 

detected in at least one other band, 
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- well fit {a^{BVR) < 0.05) by an SSP model, 

— brighter than the 90% completeness limits for each fil- 
ter in which they are detected (see Table 

1150 clusters pass these criteria, which will be used 
to study the cluster population in M51. The majority of 
the sources (from the original 3504 sources detected) that 
have been removed from our sample were eliminated be- 
cause they were fainter than one or more of our 90% com- 
pleteness limits. 

4.5. The elimination of stars on the basis of their 
energy distribution 

In i) l3.1l we argued that the expected number of very lumi- 
nous stars that might contaminate our sample is relatively 
small. We have investigated the rejection of possible stel- 
lar sources from our sample on the basis of the energy 
distribution fitting. To this purpose we created energy 
distributions of stellar atmosphere models and analysed 
them in the same way in which we analysed the energy 
distributions of clusters. 

We adopt the supergiant models calculated by Bessell 
et al. (1998), transformed to WFPC2 magnitudes by 
Romaniello et al. (2002). Bessell et al. (1998) normalized 
the magnitudes to stars with a radius of 1 Rq at a dis- 
tance of 10 pc. They pubhshed 70 models with Tgg rang- 
ing from 3500 to 50000 K for a range of gravities from 
main sequence stars to supergiants. We adopted the en- 
ergy distributions of these stellar models with the low- 
est gravity (without adding noise to their photometry) 
and then analysed the data as if the stars were clusters. 
We did this for stars that are observed in all six WFPC2 
bands {\]Y,U,B,V,R,I), and for stars observed only in 
the [/, B, V, R bands or in the B, V, R, I bands. 

We found that stars with effective temperatures higher 
than about 14,000 K all end up in the youngest age bin 
for clusters of log{t/yr) ~ 6.6. If our sample is contam- 
inated by massive stars, most of these will be blue, be- 
cause they spend most of their time on or near the main 
sequence. If there is a red phase for such massive stars, 
the crossing time between the blue and the red phase is 
very short, only about 1% of the main-sequence phase. 
However stars with masses in excess of 50 Mq are not 
expected to become red supergiants because that would 
bring them above the observed Humphreys-Davidson lu- 
minosity upper limit (Humphreys & Davidson, 1979). We 
conclude that, j/ there is a sizable number of massive O- 
stars with 21.70 < V < 22.70 mag, they will end up in the 
youngest age bin of our cluster sample. 

Fig- El shows the values of (t^{BVR) as a function of 
the temperature of the stars, for stars observed in all 
six WFPC2 bands, and for stars observed only in the 
B,V,R,I or U,B,V,R bands. Only stars below the dashed 
line pass our criterion of a^{BVR) < 0.05. We see that 
most of the stars will be rejected from our sample, except 
the coolest ones. 



25 - 




Star: T,„ ( lO' K) 

Fig. 6. The elimination of possibly contaminating bright 
stars in our sample by means of the fitting of their en- 
ergy distribution. The figure shows the values of the re- 
jection criterion a^{BVR) (Eq. 1) for stars of different 
temperatures in three different filter combinations. Notice 
that possibly contaminating stars with Toff > 17 000 K 
will be effectively eliminated by our selection criterion 
a^{BVR) < 0.05. 

We concluded before fSect. IXT)) that the expected con- 
tamination of our cluster sample of sources with V < 21.70 
mag with bright stars is very small. Here we have shown 
that our analysis of their energy distribution will elimi- 
nate a major fraction of these. The small number that 
may pass our selection criterion, are expected to end up 
in the youngest age bin. 

4.6. The effect of metallicity 

We have studied the effect of metallicity on the determi- 
nation of the age, mass, and extinction of the clusters. 
To this purpose we follow the method of de Grijs et al. 
(2003a, b) and introduce metallicity into our models as a 
free parameter. 

This technique has been applied recently to NGC 3310 
(de Grijs et al. 2003b) with good results, and to NGC 6745 
(de Grijs et al. 2003c) with somewhat less encouraging re- 
sults as two metallicities were equally well represented in 
the final sample, due to the smaller number and narrower 
wavelength coverage of available passbands, and the in- 
herent effects of the age-metallicity degeneracy. For NGC 
3310 the authors report that the near-infrared colours of 
a cluster are a good indicator of metallicity, while optical 
colours do a better job at discerning ages (see also Anders 
et al. 2004). 
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To this end we have fit all sources in M51 which pass 
our criteria (see § \AA\i with models of four metallicities, 
namely Z = 0.004, 0.008, 0.02, and 0.05. We do not find 
any significant difference in the number of sources that 
pass our selection criteria {a'^{BVR) < 0.05) between the 
different metahicity models. We find that 14%, 21%, 17%, 
and 48% of the sources in our sample are best fit (lowest 
a{BVR)^) by models of Z = 0.004, 0.008, 0.02, and 0.05 
respectively. The largest percentage of the clusters being 
fit by twice solar metallicity models is consistent with the 
metallicity estimates of HII regions in M51 by Hill et al. 
(1997). 

We note however, that many of the individual clusters 
within the same cluster associations show widely different 
metallcities. As this seems unphysical, we are forced to 
conclude that fitting on metallicity is not possible with the 
current data. We therefore choose to concentrate our anal- 
ysis on the fits using the solar metallicity models (which 
is the median metallcity found in our sample), although 
we will discuss how our results depend on the assumed 
metallicity where applicable. 
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Fig. 7. The recovery rate of our sample as a function of 
input age. We underestimate the number of clusters in 
the range of 6.9 < log(t/yr) < 7.1 by about 25% and 
overestimate the number in the preceding age bin by 
40%. 



4.7. Recovery rates 

In order to determine the recovery rate of the fitting 
method we use the artificial cluster sample generated in 
§ 14.31 In Fig. El we found that the accuracy of the fit 
depended on the magnitude of the input cluster, with 
brighter clusters being fit with higher accuracy. Because of 
the input cluster IMF, which has many more fainter clus- 
ters than brighter ones, we need to quantify this affect. 
The recovery rate is defined to be the number of clusters 
found divided by the number of input clusters for a given 
age bin. Figure shows the recovery rate as a function of 
age. This shows that we slightly over-predict the number 
of clusters in the age range of log(t/yr) G [6.62,6.87] and 
log(i/yr) G [7.37,7.62] by about 40 per cent. We underesti- 
mate the number of clusters in the age range of log(t/yr) e 
[6.87,7.12] by about 20%, and also tend to underestimate 
the number of clusters with ages greater than 10* yr. This 
is not a physical effect but due to the rapidly changing 
colours of the SSP GALEV models during this age range. 
We will use this recovery rate in S 16.21 where we study the 
cluster formation history of M51. 

5. Extinction 

Once we have derived ages, extinctions, and masses for all 
of the clusters, we can begin looking at general trends in 
the cluster population. 

Figure IHl shows the distribution of the E{B — V) val- 
ues for different age groups. The left-hand side of the fig- 
ure shows the number of clusters and the right-hand side 
shows their logarithmic values. The figure shows that at 
each age bin over half of the detected clusters have an 
extinction of E(i? — V) < 0.05 mag, with the rest show- 
ing a rapidly declining distribution up to E(i? — F) ~ 0.8 
mag. The logarithmic figures show about the same distri- 



bution of the clusters with E(i3 — V) > 0.05 mag, with 
d\og{N)/dE{B -V) ~ -2.0. This relation is independent 
of age, although the highest extinction found decreases 
strongly with ages. This is presumably due to the fading 
of clusters as they age, bringing them closer to the de- 
tection limit, so older clusters with a significant amount 
of extinction will drop below the detection limit. Another 
effect that could contribute to the observed distribution 
is if young clusters have more extinction than older clus- 
ters, which is expected as young clusters form in gas rich 
environments. 

The large numbers of clusters found with low extinc- 
tion is qualitatively what is expected. The cluster IMF of 
N{AI) ~ Af~", with a ~ 2, implies an increase in clus- 
ter numbers towards the lower mass and fainter clusters. 
The magnitudes of these low mass clusters are close to the 
detection limit, so they can only be detected if their ex- 
tinction is small. Thus our sample is complete to larger 
extinction values for brighter clusters, and for younger 
ages which are less affected by the detection limit (see 
Fig. 1101) . This effect will affect all cluster studies where 
the host galaxy is not optically thin (i.e. studies where 
one cannot see through the host galaxy at the given com- 
pleteness limits). 

In Fig. O we show the extinction as a function of age 
for the clusters in our sample. The solid points are the 
mean of all the points for each age bin and the error bars 
represent the variance about the mean. This figure shows 
that the average extinction is higher for young clusters 
(< 20 Myr) and rather abruptly drops to a constant value 
for clusters older than 20 Myr. This is qualitatively what is 
expected as clusters emerge from their parent molecular 
cloud, although as we noted above, there is a selection 
effect which causes us to detect younger clusters to higher 
extinction values. 
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Fig. 8. The histograms of the extinction for different age 
bins (given in the upper right of the panels in the left col- 
umn, in logarithmic years). The left-hand panels show the 
linear distribution {N versus Fi{B — V)) and the right-hand 
panels show the logarithmic distribution (log(A^) versus 
E{B — V)). Notice the strong peak near E{B — V) ~ mag 
and the steeply declining relation towards E(i? — ~ 0.8 
mag. 
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Fig. 9. The extinction as a function of age for the clusters 
in our sample. The solid points are the mean for each age 
bin, while the error bars are the variance about the mean. 
Note that the extinction is higher for clusters with ages 
less than 20 Myr than for the older clusters. 



6. The age distribution of the clusters 

It has been shown that the global properties of the host 
galaxy, such as an interaction or merger, can significantly 
affect the cluster formation rate (CFR) within a galaxy 



(e.g. de Grijs et al. 2003a). To search for such an effect 
within M51, we combined our cluster sample, with the de- 
tailed models of the interaction between M51 and its com- 
panion, NGC 5195, taken from the literature (SLOO). SLOO 
propose a scenario in which NGC 5195 has had two close 
passages with M51, an early encounter at approximately 
400-500 and a second at 50-100 Myr ago (the last pas- 
sage). Their double encounter model reproduces the kine- 
matic features of the system better than the single passage 
models (such as, e.g., the apparent counter-rotation of the 
southern Hl-tail), as well as the morphological "kink" in 
the northern spiral arm. 

6.1. The mass versus age diagram 

Figure irni top panel, shows the age of each cluster plotted 
against its determined mass. (We remind the reader that 
this is the present mass of the cluster.) Note the presence 
of a clear lower mass limit which increases with age. This 
is due to the evolutionary fading of clusters as they age, 
which raises the minimum mass that the clusters should 
have to be brighter than our detection limit. Hence, the 
lack of old low-mass clusters is mainly an effect of fading 
due to stellar evolution (another important effect is cluster 
disruption, which will be treated in Paper III). The solid 
lines increasing to the right show the detection limit for 
"T-F555W = 22.7 mag, i.e. the adopted magnitude limit of 
the sample, and 21.7 mag. 

The figure shows the apparent concentrations of 
the clusters in some very narrow age bins, i.e. at 
\og{t/yr)=7.2, 7.3 and near 7.8. This is not a physical 
effect but it is due to the way we select the best fitting 
models in the study of the cluster SEDs with the 3DEF- 
method (see § 14.31 and Paper I). In reality these peaks 
should spread over the neighbouring age bins. 

The bottom panel in Fig. ^] shows the relation be- 
tween the age and masses of those clusters whose sizes are 
resolved. The general properties are the same as for the 
full sample, namely the increase in the lower mass limit 
as a function of increasing age, an excess of clusters with 
young ages (< 10 Myr), and an additional excess at 60 
Myr. Therefore we conclude that we are not biasing our 
study by including unresolved sources. 

As mentioned above, certain artifacts of the age vs. 
mass distribution are likely caused by model and fitting 
effects, which do not reflect anything physical. Since this 
complicates the analysis, we attempt to correct for these 
effects using the simulations from § 14.71 To do this, we 
generate a number density grid in the age vs. mass space 
for our observed cluster sample as well as for our simu- 
lated cluster sample (cluster IMF with slope —2.0, con- 
tinuous cluster formation rate, and random extinction be- 
tween < Fi{B — V) < 1 mag with a Gaussian probability 
function that has a = 0.20 mag). The model number den- 
sity is normalized in order that the average is the same as 
the observations (i.e. we normalize to the cluster formation 
rate) . We then divide the observed cluster distribution by 



N. Bastian et al.: Stellar Clusters in M51: II. 



13 



7 F 




2 F I I I 1 

6 7 6 9 10 

log (Age/yr) 

7 t j 




2 F I I I 1 

6 7 8 9 10 

log (Age/yr) 

7 F j 



6 T 




2 F I I I 1 

6 7 8 9 10 

log (Age/yr) 

Fig. 10. Top: The present mass as a function of age for 
the 1078 clusters that pass our a criterion and which have 
a present mass greater than 1000 A/0. The solid line is the 
predicted detection limits using the GALEV models, for 
m(F439W) — 22.6 mag. The horizontal dashed lines rep- 
resent constant mass cutoffs used in Fig.^] Middle: The 
same as above, but only for those clusters with measured 
sizes greater than 2 pc. The solid line is the predicted de- 
tection Hmit for m(F555W) = 22.0 mag Bottom: The 
present mass as a function of age while leaving metallicity 
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that of the simulated distribution. This technique will re- 
move any features caused by the adopted models or the 
fitting method. 

This process and the result are shown in Fig. ^2 It 
shows the number density of the clusters in the mass ver- 
sus age diagram in a logarithmic grey-scale, normalized 
in such a way that the number densities run from light to 
dark grey. The top panel shows the density of the observed 
clusters in the mass versus age diagram. The middle panel 
shows the density of the generated cluster sample with 
a constant cluster formation rate and a power law clus- 
ter IMF. The lower panel shows the ratio of the observed 
to predicted cluster number densities in the diagram. For 
the top and middle panels, dark shading represent high 
numbers of clusters, while for the bottom dark shading 
represents an over-density of observed clusters compared 
to that expected from the models. In order to fully in- 
terpret the final result, it is important to understand the 
simulated cluster sample. Two of the strongest features 
in Fig. Illb are the increase in the number density from 
the top of the panel to the bottom as well as from left 
to right. The increase in the number of clusters at lower 
masses (top to bottom in Fig. ITTb ) is due to the cluster 
initial mass function, which gives many more low-mass 
than high-mass clusters. The increase in the number (for 
a given mass bin) of clusters from younger to older (left 
to right) is due to the logarithmic age-binning and the 
assumed constant cluster formation rate in linear time. 

Once we divide our observed distribution (Fig. Illh ^ 
by the simulated distribution (Fig. Illb ) we have the true 
age vs mass distribution which is shown in Fig. Illb . Since 
the effects of the cluster IMF and the logarithmic bin- 
ning have been removed in this comparison, the decrease 
in number density from the lower left of Fig. Illh to the 
upper right is presumably caused by the disruption of the 
clusters. Additionally we find relative excesses in the num- 
ber density at ages of ^ 60 Myr and 6 Myr. These two 
features will be discussed in more detail in § 16.21 and § 
respectively. 

6.2. The cluster formation history 

In the previous section wc found an excess of clusters at 
ages ~ 60 Myr, corresponding in time to the last passage 
of NGC 5195 and M51 in the SLOO models. In order to 
quantify this excess we will look at the age distribution 
directly (i.e. by collapsing the mass axis in Fig. Illb). As 
shown in § 16.11 for a given magnitude-limited sample, we 
are able to detect clusters to lower masses for younger ages 
than for older ages. Therefore, we will detect more younger 
than older clusters (assuming the standard cluster IMF 
which predicts many more lower mass than higher mass 
clusters), simply due to selection effects. Thus, in order to 
trace the true age distribution of the sample, the analysis 
must be restricted to masses above a certain limit, where 
the sample is complete for the age range of interest. This 
effect is shown in Fig.^lby looking at the age distribution 
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and applying the constant mass cuts shown in Fig. ^| 
The mass cuts are at 10^, 10"*, and SxlO^M©, the last 
one chosen so that we are complete for ages younger than 
1 Gyr. 

The left-hand side of Fie . II 21 shows the number of clus- 
ters detected per logarithmic age bin while the CFR (in 
number per Myr) is shown on the right-hand side. For the 
calculation of the CFR, we have corrected the age distribu- 
tion of the detected clusters for the recovery rates (Fig.[Z|) 
in order to eliminate spurious effects. The difference in 
appearance between the two sides is due to the logarith- 
mic binning of the left-hand side (e.g. each bin spans a 
larger age range than the preceding bin). Note that the 
distributions for different mass cuts appear very different. 

If we look at the sample in which we are complete for 
the past Gyr (i.e. for masses > lO*'^ Mq) we see a clear 
indication of an increase in the CFR ^ 60 Myr ago. In 
addition, we also find the peak at ^ 6 Myr, that was seen 
in Fig. lllb . This peak is presumably due to the presence of 
a large number of unbound clusters, and will be discussed 
in detail in § [T] 

The increase in the CFR at ~ 60 Myr ago for clusters 
with log(M/M0) > 4.7 occurs at the time of the proposed 
second encounter between M51 and NGC 5195. The CFR 
appears to have been 2 to 3 times higher during the burst 
period than in the period immediately preceding it. The 
CFR declines slightly after the initial rise, but remains 
significantly higher than before. In Fig. ^] we note the 
CFR rises extremely rapidly at the onset of the burst, 
and decays much more smoothly, consistent with a linear 
Gaussian distribution plotted on a logarithmic scale. 
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Fig. 11. Top (a): The present mass as a function of age 
for the 1152 clusters that pass our selection criterion. The 
shading represents the logarithmic number density, with 
the scale given on the right. Middle(b): Same as (a), 
but for the simulated cluster sample with a constant clus- 
ter formation rate (see 5 14. 31 for details). Bottom(c): The 
ratio between the observed distribution (top) and the pre- 
dicted distribution (middle). The arrows indicate two age 
bins which have an overdensity of clusters relative to the 
simulated cluster sample. 



6.2.1. Duration of the burst 

An increase in the CFR at the time of an interaction has 
been shown for other galaxies, e.g. M82, NGC 3310 and 
NGC 6745 (de Grijs et al. 2003a,b,c, respectively). While 
the duration of the burst of cluster formation in M51 is 
rather short lived, we note that this is similar to the distri- 
bution observed in NGC 3310 (de Grijs et al. 2003b) where 
the estimated duration of the increase in cluster formation 
rate was 20 Myr (it should be noted that their method 
of analysis was very similar to the one presented here). 
Presumably, the increase in the CFR in M51 was likely 
caused by the tidal interaction between M51 and NGC 
5195 60 Myr ago, supporting the two encounter model 
of SLOO. 

Tidal interactions can infiuence large portions of a 
galaxy within a relatively short time, and is the only 
known mechanism that could coordinate galaxy-wide star- 
bursts on the timescales presented here. Recently, Murray 
et al.(2004) have shown that there is an Eddington-like 
maximum luminosity for starburst galaxies. This lumi- 
nosity is a critical luminosity where there exists a bal- 
ance between gravity and radiative pressure caused by 
the young stars formed. When a starburst is tidally trig- 
gered, as is the case in M51, the (starburst) luminosity can 
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grow rapidly. When this growth happens on a shorter time 
scale than the dynamical (free fall) timescale of the star- 
bursting system, in this case giant molecular clouds, the 
luminosity grows beyond the critical luminosity and star 
formation rate stops. This could explain the narrow burst 
length observed in the cluster formation rate. Solomon et 
al.(1987) give typical values for molecular clouds in the 
Galaxy. The dynamical time scale can be estimated by: 

1/2 

idyn — Rc«/<^ — Rjs ' where a is the velocity disperion 
inside the cloud. Thus, tdyn is between 6-15 Myr for 
molecular clouds in the mass range 10^ — IO^'^^Mq as found 
in M51 (Rand & Kulkarni 1990), which are the most likely 
progenitors of star clusters in M51. This may explain the 
short duration of the burst in M51. 

6.3. Spatial Dependence of the Cluster Formation Rate 

Figure ^1 shows that the cluster formation history does 
not differ significantly between the two sides of the galaxy, 
although the west side of the galaxy seems to have ap- 
proximately twice as many clusters as the east side. This 
is not an effect of the differing detection limits between 
the two pointings, as we have restricted our analysis to 
magnitudes brighter than the 90% completeness limit of 
the shallowest observation. 

The higher cluster formation rate on the west side of 
the galaxy is also reflected in the presence of large, isolated 
and young star cluster complexes located in the western 
spiral arms, that appear to be largely absent in the eastern 
spiral arms (Bastian & Gieles 2004) . These complexes are 
typically less than ~ 6 Myr old, contain 3 x 10^ to 3 x 
1O^M0, and are associated with large concentrations of 
CO (Henry et al. 2003). 

Henry et al. (2003) studied the present day star for- 
mation rate (SFR) and its relation to the molecular cloud 
content in the spiral arms of M51, based on NICMOS 
(Paa) and BIMA-SONG (CO 1-0) observations. They find 
that both the molecular cloud content and the SFR are 
higher in the NW and W spiral arm (regions la and lb in 
their terminology) than in the SE and E arm (regions 2a 
and 2b). Thus the larger amount of clusters that we find 
in the western arm is consistent with the measured higher 
SFR on this side of the galaxy. 

7. A population of short-lived clusters 

As noted in § 16.21 the cluster formation history of M51 
shows are large increase during that last 6 Myr (see 
Fig.ll2|). This can be interpreted in two ways, either M51 
is experiencing a burst in cluster formation right now, or 
that a significant portion of the very young clusters (< 
10 Myr old) will disrupt within the next few Myr. Three 
additional observations argue for the latter scenario. First, 
a similar feature has been observed in the NGC 4038/39 
system (Whitmore 2003). Second, wc know of many loose 
open star clusters within the Galaxy that are expected 
to disperse on timescales shorter than 10 Myr (e.g. Lada 
& Lada 1991). Lastly, '-^ 10 Myr is the timescale given 
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Fig. 13. Spatial dependence of the cluster formation rate. 
We have not corrected the rates here for fitting errors. We 
note that the general trend in the formation rate is similar 
between the two sides of the galaxy. 



for NGC 5253 for clusters to dissolve and disperse their 
stars into the field (Tremonti et al. 2001). Adopting this 
scenario, we can explore what percentage of the young 
clusters will disrupt within the next few Myr. 

To do this we compare the observed "formation" rates 
in the first two bins of Fig. ^1 i.e. the bins of log(i/yr)=6.6 
to 6.93, and 6.93 to 7.26. In Fig.^]we show the percentage 
of the young clusters that are expected to disrupt before 
they age to older than 10 Myr, as a function of their mass. 
One striking feature of this figure, is that the percentage 
of the clusters disrupted is largely independent of mass 
within the uncertainties of the measurements. 

The mean value of the number ratio of clusters 
in the logarithmic age bins of (6.93,7.26) compared to 
(6.6,6.93) for clusters with a mass in the range of 3.5 < 
log(Af/M0) < 5.0 is 68± 15%. This is comparable to the 
fraction 87%) of missing clusters less than 10* years 
(observed relative to predictions based on embedded star 
clusters in GMCs) in the solar neighbourhood estimated 
by Lada & Lada (1991). This missing fraction is presum- 
ably made up of unbound clusters that have dispersed. For 
clusters to disappear from our sample due to expansion, 
their radius should extend over more than 3 pixels (with 
the exact value depending on the cluster's brightness and 
the brightness of the local background) , corresponding to 
a linear size of 48 pc. For a mean initial velocity disper- 
sion of the stars in a cluster of about 10 km s""'^, this 
corresponds to a disruption time of about 5 Myr. 

The largely mass independent destruction rate of the 
clusters is in stark contrast to other observed disruption 
patterns (Boutloukos & Lamers 2003) in which the dis- 
ruption timescale is heavily dependent on the initial mass 
of the cluster. Boutloukos & Lamers (2003) found that the 
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Fig. 12. Differences between the age distribution and the formation distribution. Left: the number of clusters found 
per age bin for three mass limits. Right: the cluster formation rate (number per Myr) as a function of age, corrected for 
the recovery errors given in Fig.El Adopting different mass limits clearly biases the interpretation of the distribution. 



disruption time depends on the cluster mass to the power 
0.6 for cluster systems in four different galaxies. Also the- 
oretical disruption studies of young cluster populations 
(e.g. Vesperini 1998, FaU & Zhang 2001, Baumgardt & 
Makino 2003) have shown a strong dependence of the dis- 
ruption timescale on the mass of the cluster. This implies 
that there are two distinct stages of disruption. The first, 
mass independent, operates on very short timescales and 
may be related to the sudden removal of gas from the 
system caused by supernovae and stellar winds (Boily & 
Kroupa 2003). This effect has been dubbed "infant mor- 
tality" (e.g. Whitmore 2004) and it corresponds to the 
rapid disruption of unbound clusters. The second stage of 
disruption, strongly mass dependent, operates on longer 
timescales and is a combination of different effects, includ- 
ing evaporation, tidal fields, and relaxation. 

Indeed, we expect that the first period of disruption 
is largely mass independent in order that the mass func- 
tions (as interpreted through the luminosity function) of 
extremely young (< 10 Myr) cluster systems such as in 
the Antennae (Whitmore et al. 1999) is similar to that of 
intermediate-age populations (300 Myr - 3 Gyr) such as in 



log (Mass/Mo 



Fig. 14. The fraction of the young (< 10 Myr old) clusters 
per mass bin that are expected to disrupt within the next 
- 10 Myr. 



NGC 7252 (Miller et al. 
et al. 2002). 



1997) and NGC 3610 (Whitmore 
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8. Cluster sizes 

Using the ishape algorithm by Larsen (1999), we have 
attempted to derive the effective radius of all clusters on 
the WF chips in M51 with an observed F439W magni- 
tude brighter than 22.6 mag. We have restricted ourselves 
to these brighter clusters, as fainter clusters do not have a 
high enough signal-to-noise ratio for unambiguous model 
fits. Briefly, the algorithm convolves a two dimensional an- 
alytical model (e.g. a King profile) with the PSF and dif- 
fusion kernel of the F555W filter for a given WFPC2 chip. 
We have modelled the PSF of the F555W filter with Tiny 
Tim, at position x — 400, y = 400 for each of the three 
wide field chips. We have tested whether the fits change 
significantly if the PSF used is from a different portion of 
the chip, and find that the difference is marginal, on the 
order of 10% in the derived radius. The algorithm then 
compares the model with the observed cluster profile us- 
ing a reduced test. For each observed cluster, multiple 
King profiles, all with concentration factor 30, of varying 
sizes were compared and the best fit model was selected. 
We then converted the measured FWHM of the profile to 
effective radius, RcS (Larsen 1999). 

Following Larsen (2004), the minimum value accepted 
for a 'resolved' cluster was a FWHM of 0.2 pixels. At the 
distance of M51, a FWHM of 0.2 pixels for a King pro- 
file with concentration factor 30 corresponds to an of 
1.2 pc. We adopt the more conservative lower limit of 2 
pc when making the subsample of resolved clusters from 
the full cluster sample. We caution that this size crite- 
rion biases the sample towards large clusters, but due to 
the resolution of the data, clusters smaller than this crite- 
rion do not have well determined radii. This results in 407 
clusters that pass the criteria given in § 14.41 have reliable 
size measurements of > 2 pc and are located more than 
3 pixels away from the nearest source. The final criterion 
was applied in order to avoid contamination by neighbor- 
ing sources in the size determination. The size criterion 
(-Reff > 2 pc) eliminates 552 clusters from the sample. Of 
the 91 sources that were detected in four or more bands, 
did not pass our a'^{BVR) criterion, and were brighter 
than V = 22.0, 84% were unresolved. This provides fur- 
ther evidence that our cluster selection eliminates the ma- 
jority of contaminating individual stars from our sample. 

Additionally, we have tested the reliability of the de- 
rived fit as a function of the magnitude of the cluster. 
To do this we generated 20 artificial clusters with ef- 
fective radii between 0.65 and 12.4 pc for magnitudes 
V = 20, 21, 22, and 23 mag (i.e., 80 clusters in total). We 
added these analytic profiles to the F555W images and 
re-fit them using the above method. Figure ITKl shows the 
results of these tests. As a function of decreasing bright- 
ness, the first clusters to deviate from the input radius are 
the large clusters. This is due to the fact that the outer 
parts of the clusters will be lost in the background noise, 
and so larger clusters will be affected first. 

The distribution of sizes for all clusters more than 3 
pixels away from the nearest source, is shown in Fig. [T?)l 
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Fig. 15. A test of the accuracy of the measurements of 
cluster radii. The derived effective radius as a function of 
both the input effective radius and the magnitude. A one- 
to-one relation is shown by the solid line. The difference 
between the recovered and the input radii is shown as the 
dashed line, while the dotted lines show the value 10% 
difference between the input and the recovered radii. 
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The dashed line in the figure is a power-law fit to the data 
of the form N{r)dr oc r^^dr, with rj — 2.2 ±0.2, where the 
fit was carried out only on clusters with R^s > 2 pc (i.e. 
the resolved sample) and Rcs < 15 pc (most larger sources 
were determined to be blends of sources upon visual in- 
spection) . We note that this fit overestimates the number 
of small clusters, a possible indication of a turnover in this 
distribution. We compare this to the size distribution of 
the Galactic globular clusters, by using the data available 
in the McMaster catalog (Harris 1996)''. Fig.llTlshows the 
size distribution (half-mass radius) of the Galactic glob- 
ular clusters, along with the best fitting power-law with 
r] = 2.4 ± 0.5. We note that the size distribution of the 
M51 clusters is different from that reported for the young 
cluster system in the merger remnant NGC 3256, which 
is characterized by = 3.4 (Ashman & Zepf 2001). In 
that same work, the authors report that the size distribu- 
tion of the Galactic globular clusters has rj = 3.4, signif- 
icantly steeper than that found here (although using the 
same catalog). This may have been caused by the erro- 
neous addition of -fl to 77 in the fitting of the size dis- 
tribution (see Elmegreen & Falgarone 1996 for a correc- 
tion of a similar point of confusion in the data of Galactic 
GMCs). We therefore report that the size distribution of 
the old Globular clusters in the Galaxy and that of the 
young M51 clusters (and possibly that of the young NGC 
3256 system) has the same form, within the errors. This 
distribution is, however, significantly different from the 
Galactic GMC size distribution, which is characterized by 
77 = 3.3 ± 0.3 (Elmegreen & Falgarone 1996), in contrast 
with what was reported by Ashman & Zepf (2001). 

The similarity between the size distributions of MW 
GCs and young star clusters suggests that disruption pro- 
cesses (which have had a much longer time to act on the 
MW GCs) have not preferentially acted upon the larger 
clusters. Thus, cluster disruption appears to be largely in- 
dependent of size. Additionally, the strong resemblance 
between size distributions of M51 and the Galactic glob- 
ular clusters lends support to the idea of a universal for- 
mation mechanism for star clusters. 
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Fig. 16. Size distribution of all clusters more than 3 pixels 
away from the nearest source. The dashed line is a power- 
law fit to the data of the form N{r)dr cx r~^dr, with 77 = 
2.2 ±0.2. The fit was carried out for all clusters with 2.0 < 
i?cff < 15.0 pc. The power-law fits the data quite well for 
sources greater than 2 pc, although it over-estimates the 
number of small clusters (unresolved clusters were given 
a size of pc) with respect to the observations. 
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9. Correlations with the cluster radius 

With a dataset such as compiled in this work, it is interest- 
ing to search for correlations among the derived parame- 
ters. In particular we are interested in relations which may 
refiect the formation mechanism of star clusters. Much em- 
phasis has been placed on the relation between the size of 
a cluster and the cluster's mass/luminosity (e.g. Zepf et 
al. 1999, Ashman & Zepf 2001). The GMCs that are pre- 
sumed to be the birth place of young clusters follow the 
relation rdoud oc M^{^^^ (Larson 1981). Young clusters, 
such as those in merging galaxies (Zepf et al. 1999) and 
in spiral galaxies (Larsen 2004) show only a slight correla- 
tion (if at all) between their radii and masses. The results 



* We used the revised February 2003 verison, 
|http://physun.physics.mcmaster.ca/Globular.html, 



Fig. 17. Size distribution (half mass radius) of the 
Galactic globular cluster population. The dashed line is 
a power-law fit to the data of the form A'^(r)dr oc r~^dr, 
with ?7 = 2.4 ± 0.5. The fit was carried out for all clusters 
with 2.0 < i?eff > 15.0 pc. 



using our sample are shown in Fig. 1181 There is no appar- 
ent direct relation between the mass and size of a cluster 
in M51. 

Recently, Larsen (2004) has shown that young clusters 
in spiral galaxies have shallower luminosity profiles than 
older clusters, presumably due to an extended halo sur- 
rounding young clusters, although no correlation between 
size and age is found. In § [S] we presented our method 
used to derive the cluster sizes, which assumed a constant 
cluster profile, unlike the study of Larsen (2004) which 
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Fig. 18. The cluster mass vs. radius for all clusters with 
i?eff > 2 pc and located more than 3 pixels away from the 
nearest neighboring source. 

also solved for the best fitting cluster profile. Figure E| 
shows the effective radii vs. ages for the resolved clusters 
in our M51 sample. No strong correlation is found, al- 
though there is a slight tendency for younger clusters to 
be more extended than their older counterparts. We also 
note that if we include the unresolved clusters (assigning 
an effective radius of 1 pc) this slight trend is removed. 

We have also attempted a multi-variable fit of the 
form: 

R^g = b X age^ * mass^. (2) 

The derived exponents, x and y, are —0.051 ± 0.027 and 
0.048±0.031, respectively, if we restrict the sample to only 
the most robust measurements, 3 < i?cff(pc) < 10 and 
sources more than 5 pixels away from the nearest source. 
As before we see a slight trend of the size decreasing with 
increasing age as well as a slight increase in size with in- 
creasing mass, though both relations are quite shallow and 
have rather large uncertainties. 

Van den Bergh et al. (1991) report that globular clus- 

1 /2 

ters in the Galaxy follow the relation Dcs oc R^i^i (with a 
large scatter) where Dcs is the diameter within which half 
the light of the cluster is contained, in projection, and -Rgai 
is the distance of the cluster to the Galactic center. They 
interpret this result as suggesting that compact clusters 
form preferentially near the centers of galaxies where the 
density of gas clouds is higher than in the outer regions of 
the galaxy. We have searched our sample in two ways to see 
if a similar effect is present. We have looked for a relation 
between the effective radius of a cluster and its galacto- 
centric distance (Fig. 120(1 as well as a trend between the 
stellar density (instead the effective radius) of each cluster 
and its galactocentric distance (Fig. 12 1|) . No clear relation 
can be seen. This may mean that the young clusters are 
forming by a different mechanism than globular clusters. 
This may simply be due to our sample coming mainly 
from the disk of M51 while the globular clusters presum- 
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Fig. 19. The age of the clusters vs. their effective radii. 
There is a slight tendency for younger clusters to be larger 
than their older counterparts, but the scatter is large. 
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Fig. 20. The effective radii of a clusters versus their dis- 
tance to the nucleus of M51. No trend is apparent in the 
data. 

ably formed outside the Galactic disk. Alternatively, the 
relation found for old globular clusters in the Galaxy may 
be due to disruption effects which have had a much longer 
time to work on globular clusters than on the young disk 
population in the M51 sample. 

10. Conclusions 

We have analysed the star cluster population of the inter- 
acting galaxy M51. By comparing the broad-band mag- 
nitudes with those of SSP models, we have derived the 
ages, extinction values and masses of 1152 clusters with 
ages between 4 Myr and 10 Gyr, and masses between 
lO^ *" — IO^Mq. The main conclusions can be summarized 
as follows: 

1. In order to examine the cluster formation history 
within a galaxy, a minimum mass cut-off must be ap- 
plied, above which the sample is complete over the 
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Fig. 21. The density of the clusters versus their distance 
to the nucleus of M51. No trend is apparent in the data. 



entire age range of interest. This avoids the bias in- 
troduced by the effects of the fading of clusters due to 
stellar evolution. 

2. The cluster formation rate increased significantly ap- 
proximately 60 Myr ago, which is coeval with the pro- 
posed second encounter between M51 and its com- 
panion NGC 5195. We interpret this as evidence 
supporting the multiple encounter model of Salo & 
Laurikainen (2000). 

3. The number of young (< 10 Myr) clusters is much 
higher than would be expected for a constant clus- 
ter formation rate over the past 20-30 Myr. By com- 
paring the formation rate during the past 10 Myr to 
that between 10-20 Myr, we find that - 70% of the 
young clusters will disrupt before they are 10 Myr old. 
This process appears to be independent of mass, which 
implies that it is a fundamentally different disruption 
mechanism from that which dominates the disruption 
of older clusters. Thus, there seem to be at least two 
distinct modes of cluster disruption. The first mecha- 
nism operates during the first ~ 10 Myr of the clus- 
ter's lifetime while the second operates on significantly 
longer timescales. We suggest that the first may be due 
to rapid gas removal from the emerging young clus- 
ter, based on the models of Boily & Kroupa (2003) 



This is similar to the fraction of unbound clusters in 
the solar neighbourhood. The second mechanism has 
been extensively modelled (e.g. Vesperini 1998, Fall & 
Zhang 2001, Baumgardt & Makino 2003) and includes 
two-body encounters within the cluster, disk shock- 
ing, and stellar evolutionary processes. This disrup- 
tion mechanism was studied for M51 by Boutloukos 
and Lamers (2003), based on the cluster study in one 
of the WFPC2 chips presented in Paper I. We will 
improve this disruption study in detail in Paper III, 
based on the new and largely extended cluster sample 
presented in this paper. 



4. The western side of M51 contains a significantly higher 
number of star clusters, although the age and mass 
distributions of the clusters present do not show any 
large discrepancies. This trend is also reflected in the 
presence of large, isolated, and young star cluster com- 
plexes located in the west spiral arm, which appear to 
be absent on the east side of the galaxy. This is also 
consistent with the estimated star formation rates of 
Henry et al. (2003). 

5. The cluster effective radius distribution can be well 
approximated by a power-law of the form N(r)dr oc 
r~'^dr, with t] = 2.2±0.2. This is remarkably similar to 
the distribution of old Milky Way Globular Clusters, 
which is characterized by 77 ~ 2.4 ± 0.5. These dis- 
tributions are significantly different from the Galactic 
giant molecular cloud (GMC) distribution {-q — 3.3). 
The agreement between the size distributions of the 
old Galactic globular clusters and the young clusters 
in M51 suggests that cluster disruption is largely in- 
dependent of cluster size and argues for a universal 
cluster formation mechanism. 

6. We have searched our dataset in order to look for corre- 
lations among the derived parameters. From this study 
we report that: 

(a) We find a strong correlation between the age of a 
cluster and its extinction, with the youngest clus- 
ters (< 10 Myr) having < Ay >~ 0.55 mag, while 
the oldest clusters (^ 1 Gyr) have < Ay 0.30 
mag. At an age of ~ 20 Myr the extinction drops 
suddenly, and then remains relatively constant for 
older ages. This is qualitatively expected as clus- 
ters emerge from the molecular cloud from which 
they formed. 

(b) A shallow relation is found between the effective 
radius of a cluster and its mass. 

(c) There is a marginal correlation between the age of a 
cluster and its effective radius, with young clusters 
being slightly larger than older ones. 

(d) Contrary to observations of globular clusters in the 
Galaxy, no correlation between the size of a cluster 
and its galactocentric radius is found. Additionally, 
no correlation between density and galactocentric 
radius is apparent in our dataset. 
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